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LOGIC FOR THE EARTH ORBITAL AEG IN THE APOLLO REAL-TIME 


RENDEZVOUS SUPPORT PROGRAM 


By Edward J. Kenyon 
SUMMARY AND INTRODUCTION 

This report presents the program logic and equations for the analytic 
ephemeris generator (AEG) that will he used in the Apollo real-time 
rendezvous support program (ARRS - ref. l) for earth orbital mission 
planning and in the Real-Time Computer Complex (RTCC). The AEG was 
initially developed in 1964 (ref. 2); the drag computation has since been 
modified. The modified drag routine and the complete flow charts are 
presented herein.' 

The development of the AEG is thoroughly described in references 2 
through 4, and the development of drag computations is presented in 
reference 3. Hence, this report only deals with changes incorporated 
into the computation of drag effects. 


Changes to the Drag Computation 


Basically, the changes to the drag computation can be broken down 
into three areas: 

1» The method by which Simpson's rule is applied has been changed. 
Previously the effects of drag were multiplied by the proper Simpson's 
rule multiplier at the end of each drag step (usually 450 seconds) until 
the total update time had been reached. Presently each drag step is 
treated independently; it is integrated as if the step were the total 
update time. This method then divides the step into two equal parts, 
which produce the three points needed for the Simpson's rule integration. 
The process is repeated until the desired time is reached. This method 
has been found to be much more accurate than the one presented in the 
original flow charts (ref. 2). 

2. In the earlier program the constants of integration, A1 and A2, 
found in the equations for the elements (e sin g) and (e cos g) in 
subroutine DRAG were computed only one time for each vector update. Now, 
to improve the calculation of the drag force, TD, these constants are 
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evaluated at the beginning of each drag step using the drag-modified 
osculating elements produced at the end of the previous integration step. 

3. The last change is in the AEG program logic. The original version 
sent the DRAG subroutine a set of osculating and mean elements at the 
beginning of the vector update, and DRAG used these elements to integrate 
out to the desired time. This method is a good one as long as the update 
is not much greater than 1 day. Longer updates should have the opportu¬ 
nity of using the AEG's complete equations so that all perturbations due 
to the earth's potential can be accurately applied to the elements. The 
change allows this to be done; the elements are now recomputed every 
2k hours by letting the AEG break the total update into 24-hour segments 
which are integrated separately by subroutine DRAG. Naturally, the 
last segment does not have to be equal to a full day. Using this improved 
method, the perturbations due to drag are added to a better set of mean 
elements, which, in turn, improves the overall prediction. 


USE OF THE AEG AND DRAG SUBROUTINES 


For the sake of brevity, the AEG and DRAG subroutines will be referred 
to collectively as "Subroutine AEG". The equations and logic have been 
programed in Fortran IV, and card decks are available for local usage. 

The argument of subroutine AEG is written (T, L 1 , Lg) and the common, 

/AEG/ (A, E, AI, G, H, AL, AN, GN, HN, R, U, K^, Kg, RD, HDP, CD, ARE, 

WHT, AIDP, RECT, SPHER, PERIOD) where 


GN. 

HN. 


a, semimajor axis, ft 
e, eccentricity 
I, inclination, rad 

g, argument of perigee, rad 

h, inertial longitude of ascending node, rad 
1, mean anomaly, rad 

n, mean motion 

g, secular rate of advance of argument of perigee, rad/sec 

h, secular rate of regression of node, rad/sec 
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HDP. 

CDi 

ARE i 

WHT. 

AIDP i 


RECT 


X 

♦ 


r, radius, ft 
u, argument of latitude 
t, time, sec 

vehicle indices, 1, 1 for the first vehicle coordinates 
only; 1, 2 for both vehicles; 2, 2 for the second vehicle 
only 

initialization control number: zero to initialize, 1 
otherwise 

drag control number: 1 if drag effects are to be computed, 
zero if not 

r, radius rate, fps 

h", secular node, rad 

Cp, drag coefficient 

A^', frontal area of vehicle, ft 2 

W L , weight of vehicle, lb 

I", mean inclination, rad 

X, Y, Z, X, Y, Z 

V, y, if;, R, X, <j> 

flight path angle 

azimuth 

longitude 

latitude 


All quantities except T, 1^, and L g are dimensioned variables, with 
two locations reserved for each. Orbital elements are input to the 
program only when they are to be initialized (K^ = 0). The drag 
constants CD., ARE., and WHT. are always input, as are K K_., L , L , 
and T. The rest are always output. i 1 2 


An example of the operations of subroutine AEG is as follows 
Initialization 

Define values (through input) for A^, E^, AI^, G^, , AL^, , 

ARE^, WHT^ (i = L^, L^), and set K.^ = 0. Call the AEG subroutine 
(K^ will be set to 1 after initialization is complete). If an error of 
input has been made (e.g. , G^ir), the message "AEG WILL NOT INITIALIZE" 
will be printed. 

Orbit Prediction 

To obtain a set of classical elements describing the vehicle's 
position and velocity at a time different from that at initialization, 
define new values for T, K gi , reset L ( and Lg, and call the AEG 

subroutine. Output will be the orbital elements and the quantities 
AN i , GN 1S HN i , FL , Ibj GD ± , Kb , HDP i , and AIDP... 
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TABLE I.- DEFINITIONS OF AEG FLOW CHART SYMBOI 


Variable 


Definition 


h 

i 

g 

h 


R L 

U 


g D- 


h 

C D 

*1 


Semimajor axis 

Eccentricity- 

Inclination 

Argument of perigee 

Longitude of ascending node 

Mean anomaly 

Mean motion 

Secular rate of change of argument of perigee 

Angular rate of change of the longitude 
of the ascending node 

Radius to either vehicle 

Argument of latitude 

Secular rate of change of argument of perigee 
due to drag 

EPOC 

Refers to Agena 
Refers to spacecraft 

Initialization control for AEG (0 for Ini¬ 
tialization, 1 for no initialization) 

Drag control for AEG 

Radius rate of change 

Mean longitude of the ascending node 
Drag coefficient 

Frontal area of vehicle 

Weight of vehicle 

Control number for reduced AEG 


Unit 

ft 

rad 

rad 

rad 

rad 

rad/sec 

rad/sec 

rad/sec 

ft 

rad 

rad/sec 

sec 


fps 

rad 

ft* 

lb 
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TABLE I.- DEFINITIONS OF AEG FLOW CHART SYMBOLS - Concluded 


Variable 


Definition 


Unit 


I" Mean inclination rad 

Subscripts 

I Initial quantity 

i Vehicle index 

j Initialization iteration index 
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TABLE II.- DEFINITION OF SUBROUTINE DRAG FLOW CHART SYMBOLS 
Variable Definition Unit 


a Semimajor axis 

e Eccentricity 

g Argument of perigee 

e sin g Eccentricity multiplied by the sine of 
argument of perigee 

e cos g Eccentricity multiplied by the cosine of 
argument of perigee 


L Argument of longitude 

U Argument of latitude 

n Mean motion number 

g Secular rate of change of argument of 

perigee 

Yg Term computed in AEG 

Constants computed in AEG 

Cg Constants computed in AEG 

C^ Constants computed in AEG 

B Sine of mean inclination 


d(e sin g) 
d(e cos g) 


Interval over which drag computed 

Mean semimajor axis 

Mean eccentricity 

Change in "g" due to drag 

Mean argument of longitude (g + l) 

Change in "e sin g" due to drag 
Change in "e cos g " due to drag 

Secular rate of change of perigee due to 
drag 


ft 


rad 


rad 

rad 

rad/sec 


rad/sec 




TABLE II.- DEFINITION OF SUBROUTINE DRAG FLOW CHART SYMBOLS - Concluded 


Variable Definition Unit 

C L Drag coefficient 

A l Frontal area of vehicle ft 2 

W L Weight of vehicle lb 

a' 1 Average mean "A" between two points ft 

dN Average rate of change of mean motion sec -1 

number due to drag 

C^q Constant computed in AEG 
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Subroutine DRAG 



Flow chart 2.- Subroutine DRAG. 
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Time step and Increment setup 



Flow chart 2.- Subroutine DRAG - Continued. 
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Mean a calculation 



Flow chart 2.- Subroutine DRAG - Continued. 
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Initial drag effects set equal to zero 



Flow chart 2.- Subroutine DRAG - Continued. 
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Q 


L = Lp+ .75y£C ? SIN C2 L D "> 
g = g+ RULE (4) 

A 1 = SINg(eSIN g) + COS g (e COS g) - C 2 ( SIN u')-C 1 (COSu 1 ) 
-y^C 51 C0S(3 u'-g) 

A 2 = COS g (e SIN g) - SIN g (e SIN g) - C 2 (COS u') + (COS u') 
-y^C 51 SIN (3 u'-g) 

C = g + (g At) 

(e SIN g) = Aj^ SIN C + A 2 COS C + C 2 SIN L + y” C 51 SIN (3 L) 

+ A 3 CCOS (g AO - 13+ RULE(5) 

(e COS g) = A x COS C - A ? SIN C + C x COS L + y" C 51 COS (3 L) 

+ A 3 CSIN CgAO]+RULE (6) 
g = T AN -1 [(e SIN g) / (e CO S g)] 

e=v4 SIN g) 2 + (e COS g) 2 
<t>= 0 



Flow chart 2.- Subroutine DRAG - Continued. 
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Calculation of eccentric anomaly 



Flow chart 2.- Subroutine DRAG - Continued. 
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Factors for variation equations 



Flow chart 2.- Subroutine DRAG - Continued. 
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Calculation of density and drag force Tp 



Flow chart 2.- Subroutine DRAG - Continued. 
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Change in orbital elements due to drag 



Flow chart 2.- Subroutine DRAG - Continued. 
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Simpson's rule section 



Flow chart 2.- Subroutine DRAG - Continued. 
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Summation of drag effects 



d a = d a + RULE (1) 
d e = d e + RULE (2) 
d (e SIN g) = d (e SIN g) + RULE (5) 
d (eCOSg) = d (e COS g) + RULE (6) 
d (gd) = d (gd) + RULE (4) 

S"=-a" + d a/2 

V 0 = \ + dij 

g Q = g 0 (V Q ) 
d v = -3/4 d a V Q /a" 

v=v 0 + d v 
9 0 = 9 0 /r > 0 
dfj= dij+ 2dr) 

V z =V o + d V 
g z = g 0 ,T ? z /’’Ig 
drj = djj + 2d^ 


\ = \ + 2d v 
a o '= a Q '+ RULE (1) 



Flow chart 2.- Subroutine DRAG - Continued. 
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Modifying and setting quantities to be returned in subroutines argument 



Flow chart 2.- Subroutine DRAG - Continued. 
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Flow chart 2.- Subroutine DRAG - Continued. 


Page 13 of 14 








22 


0 
t R = T0 + t 


V o =V' q 



Flow chart 2.- Subroutine DRAG - Concluded. 


Page 14 of 14 
































































































26 


REFERENCES 


1. Reini, W. A.: A Description of the Input to the Apollo Real-Time 

Rendezvous Support Program. MSC IN 67-FM-165, November 3, 1967. 

2. Moore, A. J.; and Reini, W. A.: Logic and Equations for Analytic 

Ephemeris Generation. MSC IN 64-FM-50, November 23, 196L. 

3. Kenyon, E. J.: Analytic Ephemeris Generation. MSC IN 68-FM-115 

May 13, 1968. 

4. Brouwer, Dirk: Solution of the Problem of Analytical Satellite 

Theory Without Drag. Astronautical Journal, Volume 6k, No. 9, 
November 1959. 





